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Abstract 

This paper answers a question asked by Ed Pegg Jr. in 2001 : "What is the max- 
imal value of min(A, B, C)/ gcd(A B, C) for A x +B y = C 2 with A, B,C > 1; 
x iVi z > 3?" Equations of this form are analyzed, showing how they map to expo- 
nential Diophantine equations with coprime bases. A search algorithm is provided 
to find the largest min/gcd value within a given equation range. The algorithm pre- 
calculates a multi-gigabyte lookup table of power residue information that is used 
to eliminate over 99% of inputs with a single array lookup and without any fur- 
ther calculations. On inputs that pass this test, the algorithm then performs further 
power residue tests, avoiding modular powering by using lookups into precalcu- 
lated tables, and avoiding division by using multiplicative inverses. This algorithm 
is used to show the largest min/gcd value for all equations with C z < 2 . 

The Tijdeman-Zagier Conjecture (T-Z Conjecture) |T||2][3][4], also known as the Beal 
Prize ProblemQ, is the claim that the following equation has no solutions: 
Tijdeman-Zagier Conjecture Equation / Beal Prize Problem Equation 
A x + B y = C z i,!/,zeZ>3 A 1 B 1 CeZ+ gcd(A,B,C) = 1 
Regarding this conjecture, in 2001, Ed Pegg Jr. asked the following question [o|: With- 
out the restriction that A, B, C be coprime, what is the maximal value of min(A, B, C) 
I gcd(A, B, C)l We define the Pegg Value to be mhi(A, B, C) / gcd(A, B, C) of a T-Z 
Conjecture Equation or Resultant Pegg Equation. 
Resultant Pegg Equation 

A x +B y = C z i,i/,zeZ>3 A,B,CeZ+ gcd(A, B, C) > 1 
Notice that in the T-Z Conjecture, the gcd(A, B,C) = 1 condition avoids an infinite 
number of trivial and uninteresting solutions. For example, adding any two integers 
A and B gives an equation, A + B = C which multiplied by A 20 B 15 C 24 yields 
{A 7 B 5 C S ) 3 + (A 5 B 4 C 6 ) i = (A i B 3 C 5 ) 5 . The Pegg Value for this construction is 
1, even if the original terms are all coprime. In contrast to this construction, any T- 
Z Conjecture counterexample will have a Pegg Value equal to min(A, B, C). And 
min(j4, B, C) will be > 1, as Preda Mihailescu showed[7) that the only perfect powers 
that differ by 1 are 2 3 and 3 2 . So in some sense, equations with a Pegg Value > 1 
are "closer" to being T-Z Conjecture counterexamples. Let us examine precisely what 
types of equations have a Pegg Value > 1 that are not T-Z Conjecture counterexamples. 



1 Original Equations and Conversion to Resultant Equa- 
tions 

Starting with a Resultant Pegg Equation with a Pegg Value > 1, then expanding and 
dividing out the common factors (so the terms to be added are coprime), one is left 
with: da x + eb v = fc z ; x,y,z > 3; a,b,c,d,e,f > 1; gcd(da,eb, fc) = 1; and 
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there exists a number TV such that Nd is an x-th power, TVe is a y-th power, and 
Nf is a z-th power. If the smallest N = 1, then the Resultant Pegg Equation is 
simply a T-Z Conjecture Counterexample that has been multiplied by a number that 
is simultaneously an x-th power, y-th power, and z-fh power. Otherwise at least one 
of d,e, f > 1, and we call the equation an Original Pegg Equation. Multiplying this 
original equation da x + eb y = fc z by TV yields its resultant equation: A x + B v = C z 
with A = {Nd) 1 /* a; B = (Ne) 1 ^ b; C = (Nf) 1 / 2 c . 

For ease of reference throughout this paper, a, b, c will refer to the bases to the 
x, y, z exponents, respectively, of the original equation. And A, B, C will refer to the 
bases to the x, y, z exponents, respectively, of the resultant equation. Similarly, d, e, f 
will refer to the original coefficients associated with the a,b,c bases, respectively. And 
D, E, F will refer to the resultant coefficients associated with the resultant A, B, C 
bases, respectively. D, E, F are the portions of the resultant bases that result from 
"spreading" the original coefficients across the equation to make each term a perfect 
power: D = (Nd) 1 / 31 = A/a; E = (TVe) 1 /*' = B/b; F = (Nf) 1 / 2 = C/c. So the 
original equation da x +eb v = fc z after multiplying by N, converts to (Da) x + (Eb) v = 
(Fc) z or equivalently A x + B y = C z .If a, b, c, d, e, / are all pairwise coprime we 
see that the Pegg Value cannot be less than min(a, b, c), because one of a, b, c will be 
a factor of min(^4, B, C) and gcd(A, B, C) will not contain any factors of a, b, c. 

Most of our original equations needing conversion to resultant form will only have a 
single coefficient > 1 and the other two coefficients will = 1 . So for ease of reference, 
when only a single coefficient > 1, that coefficient will be mapped to coefficient / 
(itself associated with base c to the z-th power) and coefficients d and e will be = 1 . 
Each of a, b, c will always have a coefficient > 1, so when a coefficient of an original 
equation is referred to as the coefficient, it specifically means the sole coefficient > 1. 

Regarding the multiplier N, consider the case that has only a single coefficient: 
±a x ± b y ± fc z = . Using the symmetrical form of the equation allows us to 
avoid repeating some arguments and also allows us to consider only positive integers 
throughout this paper. In order to convert this original equation into its resultant form, 
the equation must be multiplied by a number TV such that: TV is an x-th power, TV is a 
y-th power and Nf is a z-th power. One method of generating such an N is let TV = f q 
such that f q is an x-th power, f q is a y-th power, and f q+1 is a z-th power. We see 
that q = (mod x); q = (mod y); q + 1 = (mod z) . The Chinese Remainder 
Theorem (CRT) guarantees a solution for q with this system of congruences when 
gcd(xy, z) = 1. As an example, consider the original equation ±a x ± b v ± fc z = 
with x = y = 3; z = 5. Let TV = f q . Solving the congruences q = (mod x); q = 
(mod y); q = —1 (mod z) per the CRT gives q = 9, and multiplying the equation by 
TV = f 9 yields a conversion from the original to the resultant form: (f 3 a) 3 + (f 3 b) 3 = 
(fcf ■ 

However, an TV produced in this fashion is not necessarily the smallest number 
possible that meets those conditions. Consider the situation when the coefficient hap- 
pens to be a perfect square (such as: 1369 * 39 5 + 22505 3 = 22586 3 ) with q 2 = /, 
±<z 3 ± b 3 ± g 2 c 5 = . Instead of multiplying by the entire coefficient to the ninth 
power (g 2 ) 9 , we can instead multiply by g 3 , giving ±(ga) 3 ± (gb) 3 ± (<?c) 5 = . This 
multiplier reduction applies not only when a coefficient is a perfect power, but also to 
any perfect power that divides the coefficient. To produce the smallest TV such that TV 
is a perfect x-th power, TV is a perfect y-th power, and TV/ is a perfect z-th power, we 
need to consider the highest power of each distinct prime factor of /. 

Let p > 1 be a prime number. For an integer n > 1, the p-adic valuation v p (n) 
is defined to be the largest integer r > such that p r \ n. Consider an equation with 
d > 1; e > 1; / > 1. To convert from its original form to the resultant form, for every 
distinct prime p \ /, the entire equation must be multiplied by number p q such that q =- 
v p (d) (mod x); q =-v p (e) (mod y); q =-v p (f) (mod z) . But d, e, / are pairwise 
coprime, so any prime dividing / will not divide d or e, so the congruence reduces to: 
q = (mod x); q = (mod y); q =-v p ( f) (mod z) . When gcd(xy, z) = 1, the 
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Table 1: Minimum prime power multiplier for each prime dividing the coefficient. 
Given the exponent set {x, y, z} shown in column 1, for the original equation ±da x ± 
eb v ± fc z = 0, for each p dividing the coefficient shown in column 2, this table shows 
the minimum q for p q that the equation needs to be multiplied by to convert the original 
equation into a resultant equation. 



Exponent 
set 


Coeff associated with 
base to this exponent 


v f 

1 


(coefficient) 
2 3 4 


{4,4,3} 


3 


8 


4 1 


{5,5,3} 


3 


5 


10 | 


{3,3,4} 


4 


3 


6 9 


{5,5,4} 


4 


15 


10 5 | 


{3,3,5} 


5 


9 


3 12 6 


{4,4,5} 


5 


4 


8 12 16 


{3,4,5} 


3 


20 


40 | 


{3,4,5} 


4 


15 


30 45 


{3,4,5} 


5 


24 


48 12 36 



CRT guarantees a minimal unique solution. When gcd(xy, z) ^ 1, and gcd(c?, e, /) 
1 . the CRT guarantees that no solution exists. 

To complete the conversion of the equation from its original form to the resultant 
form, the coefficients d and e will also need to be "spread" across the resultant equation 
in a similar manner in order to make A, B, C perfect powers. 

To convert an equation from its original form da x + eb v = fc z to its resultant form, 
the entire equation must be multiplied by a number N, the smallest of which is: 

]\f — J~J pSmallost, d (p) J~J ^smallest, c (p) J~J ^smallest,/ (p) 
p\d p|e p\f 



with the functions: 

smallcst gc i(p) = smallest q where q = 
smallest^ (p) = smallest q where q = 
smallest^/ (p) — smallest q where q = 



(mod lcm(?/, z)), q =-v p (d) (mod i) 
(mod lcm(x, z)), q =-v p (e) (mod y) 
(mod lcm(a;, y)), q =-v p (f) (mod z) 



If these systems of congruences each have a solution, the CRT can guarantee a mini- 
mum product N, Table Q] shows the minimum multiplier information obtained by the 
CRT for all exponent sets with each exponent < 5. 

1.1 Non-coprime Exponents 

The CRT guarantees that an N exists to convert a given original equation when the 
moduli are pairwise coprime (in our case, when the exponents are pairwise coprime). 
But when exponents are not pairwise coprime, N can only exist when the coefficients 
of the bases of the non-pairwise coprime exponents have factors in common. But being 
that gcd(d, e, /) = 1, we know that: 

• If an original equation has two or more bases with pairwise non-coprime expo- 
nents, the equation can only be converted into a resultant form when the coeffi- 
cients of each of those bases = 1. 

• Correspondingly, the original equation cannot be converted if all three exponents 
have a factor in common (because then all coefficients = 1, and we know that at 
least one of d, e, f must be > 1). 
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1.2 x-th, y-th, z-th Power Free Coefficients 



An integer n is said to be k-free (k > 2) if for every prime p the p-adic valuation v p (ri) 
< k (that is, p k \ n). We can specify that original equation coefficients are minimized 
and their associated bases are maximized such that: d is x-th power free , e is y-th 
power free , / is z-th power free. Having an equation with v p (d) > x, v p {e) > y, 
or v p (f) > z and converting it to a resultant equation yields a final equation that is 
identical to that produced by first moving powers from the coefficients to their associ- 
ated bases such that v p (d) < x, v p (e) < y, v p (f) < z, then multiplying that equation 
by a number that is a perfect x-th power, perfect y-th power, and perfect z-th power. 
Consider this example: Exponent set {x, x, z} with gcd(a;, z) = 1; d = e = 1; and 
coefficient / is the product of non-powered prime: p\, and contains a z-th powered 
prime: p%: f = p\p\. As part of the process of converting this original equation into 
its resultant form, the entire equation must be multiplied by p\ such that p\ is an x-th 
power and p^ z is a z-th power. So: q = (mod x); q = — z (mod z) . This is 
equivalent to looking for: q = (mod x); q = (mod z) . With x and z coprime, 
the smallest q is xz. With the allowance that the entire resultant equation can always 
be further multiplied by a number, there is no reason to allow coefficients to have a 
v p (coefficient) > the exponent of the associated base. 

Putting all this information together, in order for an Original Pegg Equation to be con- 
verted to a Resultant Pegg Equation with a Pegg Value > 1, it must be of the form: 
da x + eb y = fc z x,y, z > 3 a, b, c,d,e, f > 1 

gcd(da, eb, fc) — 1 at least one of d, e, / > 1 gcd(x, y, z) = 1 
d is x-th power free; e is y-th power free; / is z-th power free; and any two bases 
that have non coprime exponents must each have an associated coefficient = 1. 
Continuing, we can also determine the types of original equations that could be con- 
verted to resultant equations with a minimum desired Pegg Value. 

1.3 Smallest Possible gcd(A, B, C) 

gcd(^4, B, C) can be no smaller than the resultant coefficient of the base(s) to the high- 
est exponent. Consider an equation with a coefficient / associated with the base to the 
z-th power. For every prime p | /, the equation must be multiplied by p q such that p q is 
an x-th power, p q is an y-th power, and p q + v p(f s > is a z-th power. The p q multiplication 
will present itself as part of D (the resultant coefficient as part of the A base) as p q ' x , 
of E (the resultant coefficient as part of the B base) as p q ^ v , and of F (the resultant 
coefficient as part of the C base) as p( q + v p(l'))/ z _ p er fJTTTJ with / > 1, we know that 
gcd(z, x) = 1, and therefore z^i. Consider the two possibilities: 

• With z < x: Because z < x, q/z always > q/x, so when q/x is an integer, 
q/z + v p (f) I z must be a greater integer. 

• With z > x: Because z > x, q/z always < q/x, so when q/x is an integer, 
q/z must be less than that integer. The sum q/z + v p (f)/z may = q/x, but it 
cannot > q/x, because in order for the sum to equal a higher integer, v p (f)/z 
must be > 1, and that is impossible (as / is z-th power free, v p (f) < z, and 
v p (f)/z < 1). 

The above logic applies to all combinations of bases and exponent sets and primes 
dividing each coefficient. An original equation must be one of the following exponent 
sets: 

• a single base has the highest exponent, and the base to the highest exponent must 
have a coefficient (such as {3,3,4}). 

• a single base has the highest exponent, and it may or may not have an associated 
coefficient (such as {3,4,5}) . 
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• two bases have the highest exponent, and each of those bases are precluded from 
having a coefficient (such as {5,5,3}). 
Regardless of the situation, the resultant coefficient of the base(s) associated with the 
the highest power will divide the resultant coefficients of the remaining base(s). 

1.4 Re-associating mm(A, B, C) 

If min(A, B,C) is not associated with a base to the highest exponent, the entire equa- 
tion can be multiplied by a number that will make mm(A, B, C) associated with a base 
to the highest exponent. For example, consider the example Original Pegg Equation 

5 3 + 427 3 = 60073 * 6 4 . 
It converts to the smallest resultant equation 

(60073 * 5) 3 + (60073 * 427) 3 = (60073 * 6) 4 
which does not have min(vl, B, C) associated with the base to the highest power. This 
equation has a Pegg Value = a = 5. Further multiplying the resultant equation by 
N = 2 lcm(x ' y ' z) = 2 12 yields 

(2 4 * 60073 * 5) 3 + (2 4 * 60073 * 427) 3 = (2 3 * 60073 * 6) 4 
which moves min(A, B, C) to C, but the Pegg Value < c as the gcd "steals" a factor 
of 2 from the c base, resulting in a Pegg Value = 3. So we do want to multiply by a 
number that will not "steal" a factor from the original base. If the resultant equation 
was instead multiplied by N = 5 12 , the final equation would be 

(5 4 * 60073 * 5) 3 + (5 4 * 60073 * 427) 3 = (5 3 * 60073 * 6) 4 
and the Pegg Value = c = 6 (at the cost of a much larger resultant equation). 

1.5 Reduced Pegg Value When Converting from Original to Resul- 
tant Form 

This "stealing" can also result when the equation is multiplied by the N to convert from 
its original form to resultant form. It can occur any time the p q as part of the resultant 
coefficient of the base to the highest power is of a smaller power than the p q as part of 
the resultant coefficient of both the other bases. As an example, Table|2]shows that this 
can occur with the {3,3,5} exponent set with v p (f) either 1 or 3. Consider the original 
{3,3,5} equation: 

(5 * 23) 3 + (2 7 ) 3 = /(3) 5 , with / = 3 * 7 * 709. 
After converting it to its resultant form: 

(/ 3 * 5 * 23) 3 + (f * 2 7 ) 3 = (/ 2 * 3) 5 
the resulting mm(A, B, C) = C, but instead of gcd(^4, B, C) being / 2 , the gcd is / 2 3, 
because it "steals" a factor of 3 from c. So the Pegg Value of the resultant equation is 
1, less than the smallest original base to the highest exponent. And there is no away 
around this "stealing" during the original to resultant equation conversion, because we 
must multiply the equation by the primes that divide the coefficient. 

1.6 Minimum Original Equation Bases for a Given Pegg Value 

Per ^1.31 gcd(A, B, C) can be no smaller than the resultant coefficient associated with 
the base(s) to the highest power. And per i ll .41 mm(A, B, C) can always be associated 
with the smallest base to the highest power. Therefore the smallest base criteria for an 
original equation to convert to a resultant equation of a given Pegg Value V is only that 
each of the base(s) to the highest exponent must be > V . Under the condition that 
the equation not be multiplied by a number to re-associate min(A, B, C) to a different 
base, min(A, B, C) / gcd(A, B, C) must immediately be > V. Consider the original 
{3,4,5} equation with x = 3, y = A, z = 5 and a coefficient of 2 associated with 
the base to the third power: ±2a 3 ± b 4 ± c 5 = 0. It converts to its resultant form: 
±(2 7 a) 3 ± (2 5 &) 4 ± (2 4 c) 5 = 0. For this equation to be able to be converted to a 
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Table 2: Power of prime in resultant coefficient after conversion from original equation 
to resultant equation. Given the exponent set {x,y, z] shown in column 1, for the 
original equation ±da x ± eb v ± fc z = 0, for each p dividing the coefficient shown in 
column 2, this table shows how the multiplier p q (with q determined from Table [TJ is 
represented in the resultant coefficients using the format [log p (D), \og p (E), log p (F)]. 



Exponent 


Coeff associated with 




i7 p (coefficient) 




set 


base to this exponent 


1 


2 3 


4 


{4,4,3} 


3 


[ 2, 2, 3] 


I 1, 1,2] 




{5,5,3} 


3 


[ 1, 1,2] 


[ 2, 2, 4] 




{3,3,4} 


4 


[1,1,1] 


[2,2,2] [3,3,31 




{5,5,4} 


4 


[3, 3,4] 


[2,2,3] [1,1,2] 




{3,3,5} 


5 


[ 3, 3,2] 


[1,1,1] [4,4,3] 


| 2, 2,2] 


{4,4,5} 


5 


[1,1,1] 


[ 2, 2, 2] [ 3, 3, 3] 


[ 4, 4, 4] 


{3,4,5} 


3 


[7,5,4] 


[14,10,8] | 




{3,4,5} 


4 


[ 5, 4, 3] 


[10,8,6] [15,12.91 




{3,4,5} 


5 


[8,6,5] 


[16,12,10] [4,3,3] 


112,9. 8] 



resultant equation with a Pegg Value > V, without further multiplying the equation to 
re-associate min(^4, B, C), the following conditions must be met: 

Da/F > V 2 7 a/2 4 > V a>V/8 
Eb/F > V 2 5 6/2 4 > V b>V/2 
Fc/F > V 2 4 c/2 4 > V c>V 

The Pegg Value may still be less than V due to gcd "stealing" when gcd(2, c) ^ 1, 
but the above minimums must be met for the equation to have a Pegg Value > V. 

2 Generating Equations with a Desired Pegg Value 

Darmon and Granville [8 j showed that with fixed coefficients and exponents, there 
can be, at most, finitely many Original Pegg Equations with unknown integers a, b, c. 
Oesterle and Masser's ABC-conjecture implies that for fixed d,e, f coefficients, even 
allowing x, y, z to vary, the total number of solutions is limited, thereby implying that 
the total number of T-Z Conjecture counterexamples is finite. 

But for our Pegg Value searching, the sizes of the coefficients on an original equa- 
tion do not have an impact on its Pegg Value - only its resultant size. As d, e, / are not 
fixed in Original Pegg Equations, the limitation proven by Darmon and Granville does 
not constrain the number of original equations. 

With W = V x+2 ~1, consider the identity {W x+2 ) x + {W x+1 ) x+l = (W X V) X+2 . 
As W and V are necessarily coprimeQ , gcd(A, B,C) = W x , and the Pegg Value 
= V whenever W X V < W x+1 . With x > 3, this is true whenever V > 1. So 
with V > 2, x > 3, the equation yields a Resultant Pegg Equation with exponents 
{x, x + 1, x + 2} and a Pegg Value of V. Therefore the answer to Pegg's question: 
"What is the maximal value of min( A, B,C)/ gcd(^4, B, C)?" is: there is no maximal 
value; we can construct an equation with any desired Pegg Value. Using the above 
method with x = 3, we can generate a resultant equation with a Pegg Value of 60000 
and a size « 2 1270 . This is far from proving that 60000 is the highest Pegg Value for 
equations < 2 1270 . Indeed, we can use a different identity on the {3,3,5} exponent 
set to generate smaller equations of a given Pegg Value. The obvious question arises: 
What is the highest Pegg Value for equations with C z < a certain size? 

1 Note that with J + K = L, any two terms being coprime forces all three to be pairwise coprime. W 
and V are coprime, because they are both coprime to 1, the third term in the equation W + 1 = V x + 2 
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3 The Highest Pegg Value in All Equations < 2 



The strategy is to guess the exponent set that will generate the highest Pegg Value 
within a given equation range. Then search the range for equations with higher and 
higher Pegg Values. Each time a new solution is found, it reduces the search space. For 
example, after finding an equation with a Pegg Value of 1000, we can exclude from 
the search any original equations that could not possibly convert to a resultant equation 
with a Pegg Value > 1000. 

After the maximum Pegg Value is found for the given exponent set within the re- 
sultant equation size range, this Pegg Value can be used to reduce (or eliminate) the 
search space of other exponent sets. A quick guess as to the type of exponent sets of 
original equations that will yield the highest Pegg Value: 

• Will have small exponents. In addition to the ABC-conjecture's expectation that 
possible solutions thin out as exponents become higher, within a given equation 
range higher exponents reduce the size of their respective base, thereby limiting 
the equation's Pegg Value. 

• Will have two identical exponents. We know that the exponent set {3,4,5} yields 
an infinite number of Pegg Value solutions, but any coefficient on this set will 
necessarily need a large multiplier to make the coefficient become a perfect 3rd, 
4th, and 5th power. This large multiplier significantly reduces the maximum 
Pegg Value within a given range. 

• Will have an exponent set where the multiplier needed to convert between the 
original equation and its resultant equation will be as small as possible. For 
example, exponent set {3,3,4} requires a square-free coefficient to be cubed to 
convert an original equation to resultant form. But exponent set {4,4,3} requires 
a square-free coefficient to be 8th powered to convert an original equation to 
resultant form. 

These conditions indicate that exponent set {3,3,4} is the best candidate exponent 
set. To gain some confidence in this assessment, we searched the abc@home database 
of 7.5 million ABC-Hits J9) looking for hits that could be formed into Original Pegg 
Equations (with each exponent < 5) that had a large Pegg Value in relationship to their 
resultant equation size. 

Recall that if non-zero positive integers A+B = C, A < B < C, and gcd(A, B, C) 
= 1, then the three are called an ABC-Triple. The radical of N is the product of the 
distinct primes dividing N. i.e., the largest square-free factor of N. The ABC-Power of 
an ABC-Triple is defined as log(C) / log(rad(ABC)). If the ABC-Power of an ABC- 
Triple is > 1, then the triple is called an ABC-Hit. 

Original Pegg Equations are not necessarily ABC-Hits. For example, 61 3 + 67 3 = 
4123 * 2 7 converts to a resultant equation with a Pegg Value of 2, yet has an ABC- 
Power of only 0.7602. But ABC-Hits are excellent candidates to have a large Pegg 
Value in relationship to their resultant equation size (because original equations with 
small radicals will tend to be those that have small coefficients). 

We define the Pegg Power to be log(Pegg Value of equation) / log(resultant equation 
size). Similar in concept to the ABC-Power, the Pegg Power provides an easy guide 
as to the Pegg Value "quality" of the equation. The highest Pegg Power found was for 
the {3,3,4} Original Pegg Equation 14 * 111 4 + 3595 3 = 3649 3 which converts to a 
resultant equation with a Pegg Value of 1 1 1 and a Pegg Power = 0.1448. This resultant 
equation had the highest Pegg Power of all ABC-Hits in the database and is shown 
in Table [4] as it has the highest Pegg Value of all {3,3,4} equations < 3558 3 . Forty- 
nine of the top fifty Pegg Powers in the database were {3,3,4} equations. The highest 
non-{3,3,4} Pegg Power in the database was thirteenth place 0.10677 for a {3,3,5} 
equation. 

This information lends support to the plan to find the highest Pegg Value in {3,3,4} 
equations < 2 100 , and use that result to reduce or eliminate the search space in other 
exponent sets. 
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3.1 Perfect Power Testing 



In order to search for the highest Pegg Value equation < 2 , we will need a perfect 
power tester for cubes, fourth powers, and fifth powers. But each power testing algo- 
rithm can be distinct. Having three different perfect power testers reduces the runtime 
inefficiency of having a general-purpose perfect power tester that checks an input for 
being a perfect y-th power against multiple y. 

Our perfect power testers work along the lines used by GMP for testing perfect 
squares iflOl . namely, before performing a rigorous but expensive y-th power test on 
an input, first verify that the input is compatible with being a y-th power modulo small 
integers. Most non-powers will quickly be eliminated, and for those inputs that pass 
the residue testing, the algorithm will then perform a rigorous perfect power test. 

For perfect cube testing, primes p = 1 (mod 3) each rule out | (p — 1) jp possible 
inputs. In addition, powers of an unused prime can assist. For cubic residues, the 
modulus 9 rules out 6/9 of possible inputs. A modulus of 27 rules out 20/27 of possible 
inputs, but being that gcd(9, 27) 7^ 1, the eliminations overlap, so it is only a minimal 
advantage to use 27 as a modulus when also using a modulus of 9. Using a somewhat 
arbitrary cut-off point, checking inputs for being compatible with being a cube modulo 
9 and the 34 primes p = 1 (mod 3) with p < 367 rules out 99.99999999999999446% 
of all inputs. 

Notice that 4 is the only exponent > 3 that needs to be checked that is not a 
prime. Checking for fourth powers against residues for primes p = 1 (mod 2) (and 
not= 1 (mod 4)) would eliminate \{p — l)/p candidates per prime, but we can do 
better by only checking against residues for primes p = 1 (mod 4) which eliminate 
|(p — l)/p candidates per prime. Checking inputs for being compatible with being a 
fourth power modulo 9, 16, 49, and the 25 primes p = 1 (mod 4) < 257 rules out 
99.99999999999999516% of all inputs. 

For fifth powers, checking inputs for being compatible with being a fifth power 
modulo 25, and the 23 primes p= 1 (mod 5) < 521 rules out 99.99999999999999571% 
of all inputs. 

Residue testing in such a manner offers another benefit. Because all our perfect 
power testing will ultimately be done on inputs of the form a x — fc z and fc z — a x , 
we do not need to compute the actual difference before performing residue tests. We 
can use modular powering, and only compute the true total when we need the rigorous 
perfect power check. For more runtime efficiency, all the modular powering operations 
can be replaced with precalculated table lookups so no actual modular powerings are 
needed. 

Our computer algorithm would instinctively seem to require expensive division 
instructions to perform residue calculations, but divisions needed for the modulus op- 
eration (% in C) can be avoided entirely by structuring the code to allow the compiler 
to optimize away actual divisions and instead use multiplicative inverses lfTTlfT2ll . In 
GCC, this requires using a hard-coded literal as a modulus (not stored in an array - 
even a constant array). 

The smallest allowable Pegg Equation exponent is 3, so using bases up to 2 64 — 1 
allows us to handle all inputs through (2 64 — l) 3 while still using the natively supported 
64-bit unsigned integer size of the host computer. This also contributes to impressively 
minimal testing time. 

As we can perform the residue tests with multiplicative inverses instead of divi- 
sions, and can perform memory lookups instead of modular powering operations, the 
time to perform a perfect y-th power test for non-y-th power inputs is effectively inde- 
pendent of the size of the inputs a, x, /, c, z, and y. 



3.2 Algorithm to Find Resultant Pegg Equations with a Minimum 
Desired Pegg Value 

This section will present an algorithm that will solve the simple case when two expo- 
nents are identical, and then conclude with a discussion of how it also works for the 
more complicated case when all exponents are different (but still have only a single 
coefficient). 

When two exponents are the same we will map them to the x and y exponents, 
allowing us to follow our convention that exponent z is associated with the base with 
the coefficient. The equation ±a x ± b v ± fc z = with x, y, z > 3; a, b, c > 1 
has three permutations when the signs are all positive. A convenient naming system is 
to title these permutations according to what b y equals: 



Equation Permutation 


Permutation Name 


a x _ j c z = b y 


ax_minus_cz 


fc z - a x = b y 


cz_minus_ax 


a x + fc z = b y 


ax_plus_cz 



When two exponents are the same, the first and last permutations listed above are 
identical, and the set of permutations to check can be reduced to ax_minus_cz and 
cz_minus_ax. 

Our algorithmic strategy is to loop through the valid / coefficients, loop through the 
valid c bases, then loop through the valid a bases then, depending on the permutation 
that is being checked, see if fc z — a x or a x — fc z yields a perfect y-th power. Consider 
the exponent set {3,3,4} with x = 3, y = 3, z = 4. There are fewer fourth powers 
within a given range than there are perfect cubes. And our perfect cube test takes about 
the same amount of time to perform as our perfect fourth-power test. So for efficiency, 
our first base loop will range through the fourth powers, and then later perform perfect 
cube testing. 

Our perfect cube test also takes about the same amount of time to perform regard- 
less of the size of the input. So regarding the two bases with identical exponents, it 
makes sense for us to loop through the range of cubes that has a smaller number of 
cubes than the other range. 

Let t = an offset > 

Let r m i n = range minimum > to search for a perfect power 

Let r max = range maximum > r m i n to search for a perfect power 
There are fewer cubes between t+r m i n and t+r max than between r m i„ and r max . For 
efficiency, we establish the convention when x = y, the a base is larger than the b base. 
With permutation type ax_minus_cz, we have a x ~ c z — b y equivalently a x = c z + b v . 
As applied to our application, there are fewer cubes between c z + b v min and c z + b y max 
than there are between b v nin and b y max . So for the permutation type ax_minus_cz, in 
comparison to looping through the b base and then checking if the sum fc z + b v is a 
perfect cube, it is more efficient to loop through the a base, then check if the difference 
a x — fcf- is a perfect cube. With type cz_minus_ax, we have fc z ~a x = b y equivalently 
a x = fc z — b y . Similarly, it will be more efficient to loop through the a base, and 
check the difference fc z — a x and see if the result is a perfect cube. 

We can extend the use of perfect power residues even further. Consider permutation 
type cz_minus_ax. For each c base, instead of looping through the entire range of 
possible a bases and seeing if fc z - a x is a perfect y-th power, we can instead loop 
through only the a bases for a given c base that are compatible with fc z - a x being 
congruent to a perfect y-th power modulo a few small integers. 

Consider again the example exponent set {3,3,4}, with x = y = 3; z = 4. and 
permutation type cz_minus_ax with / = 8; c = 103; a = 21. We know fc z — a x cannot 
be a cube because 8 * 103 4 — 21 3 (mod 7) = 2 (mod 7), which is incompatible with 
being a perfect cube (a perfect cube is = 0, 1, or 6 (mod 7)). 
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Since we know beforehand that with a fc z residue of 2 (mod 7) that an a base of 
21 is incompatible with fc z — a x being a perfect cube, there is no sense checking an a 
base of 21 in the first place. We can precalculate a large lookup array that, for a given 
fc z , shows the offsets between a bases that are compatible with the difference being a 
perfect y-th power (modulo the product of several small integers). These moduli will 
be smallest moduli from the list generated for our perfect power tester. Increasing the 
number of small moduli in this a_base_skipahead modulus product increases the size 
of the table (and table precalculation time), but decreases the number of possible a 
bases that must be checked. This table gets large quickly, because it must list the deltas 
between acceptable a bases (per the a_base_skipahead modulus) for every possible fc z 
residue. Basically, the strategy is to multiply together as many of the smallest moduli 
together to produce a single a_base_skipahead modulus product such that the lookup 
table will still fit within computer RAM (such as within a typical 4 GiB limit). Then 
any a bases that are selected for a given fc z residue will automatically be a perfect y- 
th power modulo the a_base_skipahead modulus product. So the perfect power testing 
routine will then use remaining moduli that were not used in the a_base_skipahead 
modulus product. 

In addition, sometimes we will know immediately upon checking the current fc z 
that there are no a bases that could generate a fc z — a x congruence compatible with 
the result being a perfect y-th power. For example, if fc z = 4 (mod 7), there is no a 
base such that cubing it results in fc z — a x being congruent to a cube (mod 7). In this 
case the entire a_base loop can be avoided. This lookup table is quite small because it 
is simply a list of fc z residues that can be skipped. 

The specifics of generating the a_base_elimination and a_base_skipahead tables are 
detailed in Appendix[A] The key point is that once the tables are precalculated, they 
continually eliminate over 99% of possible inputs without any further calculations. A 
single increment through the a_base_skipahead table eliminates a large swath of possi- 
ble inputs. 

The specifics of how to calculate the range of c bases, / coefficients, and a bases 
are detailed in Appendix IB! Of note is a significant reduction in search space for the 
permutation type cz_minus_ax: fc z = a x + b y . When x = y, either a x or b y must be 
> fc z /2, and the other must be < fc z /2. This reduces the range to check from the 
standard range: a x = a x mm ..fc z - b v min [when x ^ y] to a x = fc z /2..fc z - b v mm 
[when x = y]. With this reduction, it takes much less time to check the cz_minus_ax 
permutation than the ax_minus_cz permutation. 

Notice that the gcd condition on line|7]is checked after first verifying that the appro- 
priate difference is a perfect power. If the gcd check was performed before checking for 
a perfect power, it would rule out 29% of possibilities (approximately (1 — (6/ir 2 )) 
|[T3l ) but even though gcd is an efficient algorithm, we can test for perfect cubes even 
more quickly than calculating gcd when most inputs are not perfect cubes. And being 
that the difference is rarely a perfect cube, it is more efficient to test the gcd condition 
after first verifying the perfect y-th power condition. 

A note on line [8j "Calculate maximum Pegg Value with this original equation 
within given resultant equation range" 

The first step is to multiply the entire equation by smallest N to produce a valid 
Resultant Pegg Equation. At this point, the Pegg Value will be > the desired Pegg 
Value unless either: 

• min(A, B, C) is not associated with the base to the highest exponent. In this 
case, the Pegg Value may be able to be increased by multiplying the entire equa- 
tion by a number that moves mm(A, B, C) to a different base. Of course, the 
new resulting equation still must fit within the desired equation range. See fll.41 

• min(^4, B, C) is associated with a base to the highest exponent, and the gcd 
"steals" a factor from the base associated with the original coefficient. In this 
case, nothing can increase the Pegg Value. See fll.51 
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Algorithm 1 Within the specified equation range, find an equation with a minimum 
desired Pegg Value that is based on an Original Pegg Equation with a single coefficient 



1: for permutation in ax_minus_cz to cz_minus_ax do 

2: for c base in all possible inside resultant eq range with this permutation do 
3: for / coeff in all possible inside resultant eq range with this permutation and 
c base do 

4: if a_base_elimination table does not rule out all a bases for this / and c z 

then 

5: for a base in all possible inside resultant eq range with this permutation 

and c base and / coeff [use the a_base_skipahead lookup table to only 
check the a bases that are compatible with the result being congruent to 
a y-th power modulo a_base_skipahead modulus] do 
6: if (permutation = ax_minus_cz and a x — fc z = perfect y-th power) 

or (permutation = cz_minus_ax and fc z — a x = perfect y-th power) 
then 

7: if gcd(a, c) = 1 then 

8: Calculate maximum Pegg Value with this original eq within 

given resultant eq range 
9: if the maximum Pegg Value > minimum desired Pegg Value 

then 

10: return Pegg Value,Resultant Pegg Equation 

ii: end if 

12: end if 

13: end if 

14: end for 

15: end if 

16: end for 

17: end for 

18: end for 

19: return entire range searched - none exist 



This algorithm can also be used for exponent sets where all exponents are different by 
reordering the input parameters. The standard input parameters are {x, y, z} where x 
and y are exponents associated with bases with no coefficients. 

z = exponent associated with the base with coefficient (the algorithm will loop 
through these c bases) 

x = exponent associated with the base with the smaller number of elements (the 
algorithm will loop through these a bases) 

y = exponent associated with the base with the larger number of elements (these 
are detected by performing perfect y-th power tests on fc z — a x or a x — fc z ) 

When x = y, the permutation ax_plus_cz was ignored, as it was identical to az_minus_cz. 
But when x ^ y, the unprogrammed permutation ax_plus_cz can be checked by re- 
ordering the inputs to the algorithm. Consider the exponent set {3,4,5} where the 
coefficient is associated with the base to the third power. We follow our convention 
that exponent z is associated with the base with the coefficient, so z = 3. 

The following table shows how we check all three permutation for z = 3 with our 
current algorithm which handles the two permutation types ax_minus_cz and cz_minus_ax: 



Equation Permutation 


Handled with input parameter ordering {x, y, z} 


a x - fc z = b v 


ax_minus_cz with exponent set {5,4,3} 


fc z - a x = bv 


cz_minus_ax with exponent set {4,5,3} or {5,4,3} 


a x + f c z = by 


ax_minus_cz with exponent set {4,5,3} 
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Table 3: All combinations of T-Z Conjecture counterexample exponent sets with each 
exponent < 5, along with a possible reference to a proof showing that none exist 



exponent 


reference to proof that exponent set has 


set 


no T-Z Conjecture counterexamples 


{n,n,n} 


Wiles and Taylor. Fermat's last Theorem. 118111911 




although the specific cases with 3 < n < 5 were all solved earlier 


{3,3,3} 


Euler (18th century) 


{4,4,4} 


Fermat (17th century) 


{5,5,5} 


Dirichlet and Legendre (19th century) 


{3,3,4} 


Bruin 1 1 1 


{3,3,5} 


Bruin 1 1 1 


{4,4,3} 


Lucas (19th century) 


{4,4,5} 


Bruin |20| showed complete list of {2,4,5} 


{5,5,3} 


Poonen Uli 


{5,5,4} 


Poonen [21 1 ruled out {5,5,2} 


{3,4,5} 


has not been ruled out 



Notice that the second possibility can be checked with exponent set {4,5,3} or 
exponent set {5,4,3}. Our ax_minus_cz and cz_minus_ax search algorithms run more 
efficiently with x > y, so it will be most efficient to check exponent set {5,4,3}. 

For the third possibility, it is less efficient at runtime to check ax_minus_cz with 
exponent set {4,5,3} than using a new ax_plus_cz routine and checking exponent set 
{5,4,3}, but the inefficiency is not great enough to overcome the overhead of writing 
a new routine. This searching does not take very long because there are far fewer 
combinations to check when all three exponents are different. 

Using the above algorithm, it is easy to determine that the smallest {3,3,4} equation 
with a Pegg Value > 1 is 207 3 + 126 4 = 639 3 which has Pegg Value of 14. Searching 
{3,3,4} equations up to 2 100 for higher and higher Pegg Values yields the information 
in Tableland FigureQ] The last row in Table|4]shows a resultant equation < 2 100 with 
a Pegg Value of 63742. The next step is to prove the hunch that this is the highest Pegg 
Value in all equations < 2 100 . 

Lemma 1 There are no T-Z Conjecture counterexamples < 2 100 with a Pegg Value 
> 63742 

Any equation with a higher Pegg Value must have all bases > 63742. This pre- 
cludes any base having an exponent of 7 or higher, as 63743 7 > 2 100 . 

A good summary of exponent sets for which it has been proven that no T-Z Con- 
jecture counterexamples exist is provided in both H141 and fTBI . Of all exponent sets 
with each exponent < 5, all but one have been ruled out as having any T-Z Conjecture 
counterexamples. This information is summarized in Table[3] Note that exponent 6 is 
covered under exponent 3, as a 6 = (a 2 ) 3 , so for T-Z Conjecture counterexamples we 
only need to consider exponent sets with each exponent < 5. 

The only T-Z Conjecture counterexample exponent set listed above that has not 
been ruled out is {3,4,5}. Notice that {3,4,5} can be rewritten as a special form of 
{2,3,5} for which the complete parameterization was shown by Johnny Edwards |[T6l 
and conveniently available on Dario Alpern's web-site [17]. Using these parameteriza- 
tions, checking all solutions for 

±a 2 ± b 3 ± c 5 — with a,b,c> 1, gcd(a, 6, c) = 1 
shows that for equations < 2 800 , in no case is a, the base to be squared, itself a perfect 
square. So there are no {3,4,5} T-Z Conjecture counterexamples with an equation 
size < 2 800 , thereby proving that there are no T-Z Conjecture counterexamples of any 
exponent set under an equation size < 2 100 that have a Pegg Value > 63742. 
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Table 4: Smallest {3,3,4} resultant equation with a Pegg Value higher than that shown 
in the preceding row (Row 1 shows the smallest {3,3,4} equation with a Pegg Value 
> 1). 



log 2 








(resultant 


Pegg 


Pegg 




equation) 


Value 


Power 


Original Equation 


27.96 


14 


0.1362 


23 a + 9 * 14 4 = 71 a 


33.81 


21 


0.1299 


13 * 21 4 + 163 ;1 = 190 a 


43.80 


43 


0.1239 


23 * 43 4 + 1056 3 = 1079 3 


46.92 


111 


0.1448 


14* 111 4 + 3595 3 = 3649 3 


56.75 


133 


0.1243 


1157 a + 139* 133 4 = 3558 a 


57.82 


183 


0.1300 


1966 3 + 121 * 183 4 = 5233 a 


60.68 


194 


0.1252 


126 * 194 4 + 9071 a = 9743 a 


66.96 


201 


0.1143 


5906 a + 8809 a = 545 * 201 4 


66.98 


365 


0.1271 


10973 a + 15902 a = 301 * 365 4 


69.24 


399 


0.1248 


12146 a + 391 * 399 4 = 22703 a 


72.75 


455 


0.1214 


513 * 455 4 + 33247 a = 38872 a 


73.74 


1482 


0.1429 


1609 a + 239 * 1482 4 = 104857 a 


73.81 


1638 


0.1447 


97103 a + 193 * 1638 4 = 132095 3 


74.25 


2994 


0.1555 


104 * 2994 4 + 226199 3 = 271127 a 


90.12 


3858 


0.1322 


25031 a + 1570 * 3858 4 = 70327T 5 


90.16 


5838 


0.1388 


729217 a + 971 * 5838 4 = 1148689 3 


90.82 


11598 


0.1487 


341 * 11598 4 + 3662591^ = 3809903 a 


92.75 


49476 


0.1681 


7771657 a + 8824055 a = 193 * 49476 4 


99.91 


63742 


0.1597 


2192137 a + 20440855 3 = 518 * 63742 4 




1 L- ' ' ' ' ' ' ' 1 

30 40 50 60 70 80 90 100 

log 2 of resultant equation size 

Figure 1: Highest {3,3,4} Pegg Values for a given resultant equation size (information 
from Table 3|i 
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Lemma 2 There are no Resultant Pegg equations < 2 100 with a Pegg Value > 63742 
with an exponent of 6 

In an equation < 2 100 , the maximum possible base to a sixth power is |_2 100//6 J = 
104031. The only reason to use exponent 6 instead of 3 on a particular base would 
be if one wanted to force that particular base to be smaller so that it would be se- 
lected as min(A, B, C). When mm(A,B,C) is associated with an exponent of 6, 
mm(A,B,C)/ gcd{A,B,C) can be no higher than [104031/2J = 52015, because 
gcd(^4, B, C) must be > 2. As 52015 < 63743, there is no reason to consider equa- 
tions with an exponent of 6. 

Lemma 3 There are no Resultant Pegg equations < 2 100 with a Pegg Value > 63742 
with two exponents the same 

The Original Pegg Equation condition that gcd(a;, y, z) = 1 precludes z = x, and 
as all exponents to consider are < 6, the following sets must be checked: {4,4,3}, 
{5,5,3}, {3,3,4}, {5,5,4}. {3,3,5}. {4,4,5}. Regarding the search for a {3,3,4} equa- 
tion with a coefficient of 2, Henri Cohen l22l showed that the equation: 

±a 3 ±b 3 ± 2c 2 = 0, with a, 6, c > 1, gcd(a, b, c) = 1 
can be parameterized by one of four parameterizations. Using these parameterizations, 
checking all solutions with 2c 2 < 2 240 shows that in no case is c, the base to be 
squared, itself a perfect square. So there are no {3,3,4} original equations < 2 240 
with a coefficient of 2, thereby eliminating any need for our algorithm to consider a 
coefficient of 2 for exponent {3,3,4} in a search for any resultant equations < 2 100 . 

The remaining coefficients for {3,3,4} and the other listed exponent sets were 
checked with Algorithm Q] and in no case was an Original Pegg Equation found that 
could be converted to a Resultant Pegg Equation < 2 100 with a Pegg Value > 63742. 

Lemma 4 There are no Resultant Pegg Equations < 2 100 with a Pegg Value > 63742 
with no two exponents the same 

As all exponents to consider are < 6, the only the set to check is {3,4,5}. All 
the other exponent sets to check had two identical exponents, so we mapped the non- 
identical exponent to z according to our convention that the sole base that had a coef- 
ficient > 1 would be mapped to the c base. For exponent set {3,4,5} all bases could 
have a coefficient > 1, so it is not obvious how to map the x, y, z exponents. To ease 
confusion, for now we will simply map: x = 3; y = 4; z = 5. 

First to note is that there is not room to multiply the entire equation by a number in 
order to increase its Pegg Value by moving min(v4, B, C) from being associated with a 
smaller exponent to a larger exponent: Per fll.31 the smallest possible gcd(j4, B, C) = 
F, so the largest Pegg Value for a {3,4,5} equation is min(A,B,C)/F. The high- 
est value for C is L2 10 °/ 5 J = 1048576. So min(A, B, C) must be < 1048576. As 
the Pegg Value is limited to mm(A, B, C)/F, the Pegg Value < [1048576/FJ. For 
the Pegg Value to be > 63743, F must be < L1048576/63743J < 16. The smallest 
number by which the entire equation could be multiplied is 2 lcul ^ x ' v ' z ^ = 2 60 , which 
would contribute to a resultant coefficient F of 2 60 / 5 = 40 96, which is greater than the 
maximum of 16. 

Second to note is that {3,4,5} can have only single coefficient. There must be at 
least one original coefficient > 1, therefore the resulting D, E, F will all be > 1, and 
we also know that F < 16. We must answer the question, "Under what conditions can 
an original {3,4,5} equation result in 1 < F < 16?" 

Our search is limited to F < 16 < 2 4 , so v p (F) < 4, for every p | F. Looking at 
Table|2]for options where a {3,4,5} equation has F < 4th power, we have: 

{3,4,5}, with v p (d) = 1, F = p i {3,4,5}, with v p (e) = 1, F = p 3 

{3,4,5}, with v p (f) - 3, F = p 3 
There cannot be coefficients on more than one base (because there is only room for a 
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maximum resultant coefficient of 2 4 for F and the smallest combination of two coprime 
coefficients would yield a resultant coefficient of F = p\p\ which would be > 2 4 , the 
maximum value of F). And p can only = 2, as the smallest possible F is p 3 and cubing 
the next higher base yields 3 3 which is greater than the maximum F of 16. The search 
for {3,4,5} original equations that convert to resultant equations < 2 100 with a Pegg 
Value > 63742 is limited to: 

{3,4,5} with a coefficient of 2 associated with the base to the 3rd power. 

{3,4,5} with a coefficient of 2 associated with the base to the 4th power. 

{3,4,5} with a coefficient of 8 associated with the base to the 5th power. 
As the equation cannot be multiplied by a number to re-associate mm(A, B, C), £ 11.61 
allows us to calculate the original base minimums for each base (and shows the first 
{3,4,5} case worked as an example). The three coefficient scenarios were all checked 
with Algorithm Q] and in no case was an Original Pegg Equation found that could be 
converted to a Resultant Pegg Equation < 2 100 with a Pegg Value > 63742. 

Conclusion 

The last row of Table|4]shows a {3,3,4} Original Pegg Equation that can be converted 
to the equation 1135526966 3 + 10588362890 3 = 33018356 4 that has a Pegg Value of 
63742 and a C z < 2 100 . Lemmas[T]through|4]show that there are no equations < 2 100 
of any exponent set that have a higher Pegg Value. 

Searching and proof verifying took place over several weeks on three AMD 64-bit 
computers; from two to four cores - 2.5 Ghz to 3.0 Ghz; 4 GiB to 8 GiB of memory; 
running NetBSD 4.0, 4.0.1 or 5.0 03 with the operating system included GCC 4.1.2 
or 4.1.3 and using GMP 4.2.3 or 4.3.1 ESI . 

The highest Pegg Power of the equations in Table [4] is 0.1681 . The Pegg Power 
limit for {3,3,4} equations and {4,4,3} equations is limy^oo log(V)/ log((2T^) 4 ) = 
0.25. Any T-Z Conjecture counterexample or other Pegg Equation exponent set will 
have at least one exponent > 5, which limits the Pegg Power to 0.2. So the highest 
Pegg Power for all equations lies somewhere within the range 0. 168 1 and 0.25. Rather 
than simply looking to increase the lower limit by further searching, a better technique 
to reduce this range could help shed light on reducing or eliminating possible T-Z 
Conjecture counterexamples. What is the highest Pegg Power? 
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5 Electronic Availability 

The Pegg Value searching code, written in C with calls to GMP, is released under the 
GNU General Public License (GPL) and is available at 
http://www.arthurtownsend.com/peggvalue.htm 

A Generating the a_base_skipahead and a_base_elimination 
Lookup Tables 

This appendix details the suggested a_base_elimination moduli and a_base_skipahead 
moduli for all exponent sets with each exponent < 5. The use of these tables was 
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outlined in ^3.2\ 

The a_base_elimination table is a list of fc z residues for which there are no a bases 
that can generate a fc z and a combination that is a y-th power. For these fc z , the entire 
a_base loop can be avoided. 

The a_base_skipahead table is a large lookup array showing the offsets between pos- 
sible a bases for a given fc z (offsets between a bases that are not incompatible with 
the fc z and a x combination being a perfect y-th power modulo the a_base_skipahead 
modulus). Any a bases that are selected will automatically generate an fc z and a x 
combination that is congruent to a perfect y-th power (modulo a_base_skipahead mod- 
ulus), so the perfect power testing routine will then use y-th power testing moduli that 
were not used in the a_base_skipahead modulus product. 

The tables based on these moduli were used to generate the information in Table [4] 
and also used for verifying Lemmas [3] and 2J. 

The bulk of the time of the algorithm will be spent on the coefficients that pro- 
duce the largest number of combinations of c bases and a bases within a given resul- 
tant equation range. For {3,3,4}, these will be the smallest coefficients. The {3,3,4} 
a_base_skipahead table for all possible coefficients is 2.6 GiB for the modulus 13 * 19 
* 31 * 37. The single coefficient "3" provides so many c base and a base combina- 
tions within a given range that it is worthwhile to generate a {3,3,4} single-coefficient 
a_base_skipahead table just for the coefficient 3. For a given a_base_skipahead mod- 
ulus, having an a_base_skipahead table for only a single coefficient produces a much 
smaller table than a table that can handle all coefficients. And because the resulting 
table is much smaller, we can use a larger modulus 7*13*19*31* 37 and keep 
the table size less than our 4 GiB limit. Using the additional 7 as part of the modulus 
rules out 2/3* (7— l)/7 « 57% more a bases than the original a_base_skipahead ta- 
ble. Note that these savings do not result in 57% fewer a bases being checked, because 
some of those now rejected would already have been rejected by the a_base_elimination 
modulus which also includes the modulus of 7. But the benefits of a custom individual 
coefficient a_base_skipahead modulus are sufficiently great that they can be used for a 
few coefficients on the {3,3,4} exponent set (as shown in Tabled column 4). It is not 
worthwhile using a custom a_base_skipahead table for every coefficient on the {3,3,4} 
exponent set, because the higher the coefficient the fewer combinations it produces, and 
the benefit of the 57% a base reduction must be balanced against the single-coefficient 
table precalculation time. 

Some notes on the tables: 

Regarding the coefficients tested individually (listed in Table|6] Column 4): 

• TableQ]shows that regardless of the powers of the prime factors that divide /, for 
exponent set {3,3,4} the smallest N possible is always N = f 3 . So for {3,3,4}, 
the smaller the coefficient, the more combinations there are to be checked. 

• Table[T]shows that the multiplier to convert a {4,4,3} Original Pegg Equation to 
its resultant equation is smallest when v p (f) = 2, so for {4,4,3}, the smaller and 
the "more square" the coefficient, the more combinations there are to be checked. 

Exponent set {3,3,4} shows the coefficient 7 being included in the list of coeffi- 
cients checked individually. For coefficient 7, the listed individual a_base_skipahead 
modulus 7*13*19*31* 37 offers no reduction in the a base search space over the 
smaller 13 * 19 * 31 * 37. So for coefficient 7, either the standard a_base_skipahead 
table or single-coefficient a_base_skipahead table could be used with identical effect, 
or searching efficiency could be improved by using a modulus of 9 * 13* 19*31 * 37. 

Exponent set {3,3,5} does not have coefficients checked individually because after 
increasing the standard a_base_skipahead modulus, even for use with just a single coef- 
ficient, the resulting {3,3,5} single-coefficient a_base_skipahead table is 16 GiB (which 
exceeds the 4 GiB limit). Other exponent sets do not take long to check, so it is not 
worth the overhead to precalculate single-coefficient a_base_skipahead tables for them. 

For a given exponent set, the ax_minus_cz and cz_minus_ax a_base_elimination 
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Table 5: Suggested a_base_elimination moduli and standard a_base_skipahead moduli. 

For each exponent set {x, y, z} shown in column 1 and permutation type shown in column 2, column 3 shows the a_base_elimination modu- 
lus to generate the list of fc z residues that result in no possible a x being able to produce a fc z and a x combination that is perfect y-th power. 

Column 4 shows the standard a_base_skipahead modulus. This modulus results in an a_base_skipahead table less than 4 GiB that can handle 
all coefficients. 

For all possible c bases and for all / coefficients < 100000, columns 5, 6, and 7 show how many a bases are eliminated by the listed lookup table 
as being able to produce an fc z and a x combination that could be perfect y-th power. Column 5 shows what percentage would be eliminated solely 
based on the a_base_elimination modulus. Column 6 shows what percentage would be eliminated solely based on the a_base_skipahead modulus. 
Column 7 shows what percentage were eliminated based on both tables. 













% of inputs 


% of inputs 


% of inputs 










standard 


eliminated by 


eliminated by 


eliminated by 


exponent 


permutation 




a_base_elimination 


a_base_skipahead 


elimination 


skipahead 


combination 


set 


type 




modulus 


modulus 


modulus 


modulus 


of moduli 


{3,3,4} 


both 




7*3 2 


13* 19*31*37 


47.149 


97.596 


98.729 


{3,3,5} 


both 




7*3^ 


13* 19*31*37 


46.956 


97.596 


98.725 


{4,4,3} 


ax_minus_cz 




5 * 2 4 * 17 


3^*13*29*37 


77.398 


98.063 


99.562 


{4,4,3} 


cz_minus_ax 


5*13*2 4 


* 17 * 3 a * 29 *7 2 * 11^*43^ 


3 2 * 13* 29* 37 


91.228 


99.160 


99.853 


{4,4,5} 


ax_minus_cz 




5 * 2 4 * 17 


3^*13*29*37 


69.265 


98.022 


99.392 


{4,4,5} 


cz_minus_ax 


5*13*2 4 


* 17 * 3 a * 29 *7 2 * 11^*43^ 


37*41 *53 


87.963 


98.042 


99.764 


{5,5,3} 


both 




11*5^*31*41*61 


61 * 71 * 101 


87.392 


98.915 


99.830 


{5,5,4} 


both 




11*5^*31*41*61 


41 * 61 * 71 


87.332 


98.801 


99.767 


{4,3,5} 


both 




13 


7* 3^ * 19 * 31 


7.101 


99.475 


99.512 


{5,3,4} 


both 




31 


7* 3^ * 13 * 19 


49.951 


99.773 


99.886 


{3,4,5} 


both 




13 


5*3^*2 4 *17*29 


7.101 


99.740 


99.758 


{5,4,3} 


both 




11 *41 


5 * 3^ * 13 * 2 4 


17.002 


99.638 


99.699 


{3,5,4} 


both 




31 


11*5^*41*61 


49.951 


99.993 


99.997 


{4,5,3} 


both 




11 *41 


5 2 * 31 * 61 


17.002 


99.865 


99.888 



Table 6: Suggested single-coefficient a_base_skipahead moduli. For each exponent 
set {x, y, z} shown in column 1 and permutation type shown in column 2, column 3 
shows the single-coefficient a_base_skipahead modulus. This modulus results in an 
a_base_skipahead table less than 4 GiB that can handle just a single coefficient. For the 
two exponent sets that take the longest to search, it is worthwhile to check the coeffi- 
cients that result in the most combinations with their own individual a_base_skipahead 
table. Column 4 lists the coefficients that were checked individually to generate the 
information in Table [4] and for verifying Lemmas [3] and [4] The list of coefficients in 
column 4 is sorted by decreasing number of combinations the coefficient produces. 







single-coefficient 


coefficients 


exponent 


permutation 


a_base_skipahead 


checked 


set 


type 


modulus 


individually 


{3,3,4} 


both 


7 * 13 * 19 * 31 * 37 


3,4,5,6,7,8,9,10,11,12,13 


{4,4,3} 


both 


5*3^*13*29*37 


4,2,9,25,3,36,49,18,100 



moduli are different only when two exponents are even, which when each exponent 
< 5 means when two exponents = 4, which means {4,4,3} and {4,4,5}. 

When based upon identical a_base_skipahead moduli, the generated a_base_skipahead 
tables for ax_minus_cz and cz_minus_ax permutation types for each exponent set are 
identical except when y is even, which when each exponent < 5 means when y = 4, 
which means {4,4,3}, {4,4,5}, {3,4,5}, {5,4,3}. 

The a_base_skipahead table is overwhelmingly size-sensitive to the original mod- 
uli size, so for non-prime moduli typically the smallest helpful prime power is used 
(such as 9 for {3,3,4} and {4,4,3}), even though a higher power (such as 27) elimi- 
nates slightly more candidates. But the a_base_elimination table is not under such strict 
size restrictions, so it can easily use the higher power if the higher power assists. For 
{3,3,4} the a_base_elimination modulus 27 offers no advantage over 9, but for {4,4,3} 
cz_minus_ax, the a_base .elimination modulus 27 does eliminate a few more candidates. 

The general plan is to use whatever small moduli as part of the a_base_elimination 
modulus that help eliminate a bases, then use small moduli coprime to those as part of 
the a_base_skipahead modulus. This combination of coprime moduli typically elimi- 
nates the most candidates. But some exceptions are noted: Exponent set {4,4,3} per- 
mutation type cz_minus_ax uses the modulus 29 as part of both the a_base_elimination 
modulus and the a_base_skipahead modulus. Instead of using 29, if the a_base_skipahead 
modulus was to use the next smallest unused perfect y-th power modulus (41) the re- 
sulting a_base_skipahead table would be > 4 GiB. Keeping the modulus of 29 in the 
a_base_elimination modulus provides for quicker elimination for those fc z residues 
that do not have any valid a base combinations because the a_base_elimination ta- 
ble is more likely to be stored in cache. And keeping the modulus of 29 in the 
a_base_skipahead table allows the a_base_skipahead modulus to reduce the number of a 
bases for those fc z residues that do have valid a bases. Similarly, {4,4,3} permutation 
type cz_minus_ax uses 3 2 as part of the a_base_skipahead modulus and, as mentioned 
in the previous item, uses 3 3 as part of the a_base_elimination modulus. 

The a_base_elimination modulus may not exclude as many possibilities as possi- 
ble. For example, for {4,4,3} permutation type cz_minus_ax, the a_base_elimination 
modulus contains a modulus of 16, even though 1024 would rule out a few more can- 
didates. The small number of candidates that would be excluded by a modulus of 
1024 that are not excluded by the modulus of 16 are not worth the overhead of the 
increased a_base_elimination table size. Instead, these few candidates that make it past 
the a_base_elimination check will be immediately rejected by the a_base_skipahead ta- 
ble. 

Exponent set {4,4,5} permutation type cz_minus_ax could use the same a_base_skipahead 
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modulus as {4,4,5} permutation type ax_minus_cz, but the {4,4,5} permutation type 
cz_minus_ax a_base_skipahead modulus is modified to be the product of moduli that 
are coprime to the exponent set's a_base_elimination modulus. 

Exponent set {5,5,4} could use the same 61 * 71 * 101 a_base_skipahead modulus 
as exponent set {5,5,3}, but set {5,5,4} takes such a short time to check all possible 
equations that it is not worth the extra table precalculation time. 

All the a_base_skipahead table offsets fit within a 16-bit unsigned table entry (all 
offsets < 65535), other than some of the {4,4,3} single-coefficient a_base_skipahead 
tables. For some {4,4,3 } single-coefficient tables (such as for coefficient 3, permutation 
cz_minus_ax), the a_base_skipahead modulus eliminates so many a base candidates that 
sometimes the delta between valid a bases for a given fc z residue is > 65535, requiring 
the entire table to use 4-byte fields instead of the 2-byte fields required for all the other 
tables. 

As mentioned in £13.21 the a_base_elimination lookup table is quite small because it 
is simply a list of fc z residues that can be skipped. The list could be stored through a 
variety of methods (such as a hash table or binary search tree), although storing it as 
a simple binary array provides for quicker access. A binary array requires a number 
of elements equal to the modulus product. With a sufficient number of primes in the 
modulus product even this table can become large, in which case the modulus product 
can be split into different moduli (into multiple lookup tables), with each table able 
to reduce possibilities (unlike the a_base_skipahead table which must be just a single 
table). Even though Table [5] shows only a single a_base_elimination modulus for each 
exponent set, for the purposes of code implementation, the moduli for exponent sets 
{4,4,3 } and {4,4,5} permutation cz_minus_ax were each split into three to provide three 
distinct a_base_elimination lookup tables. And the moduli for {5,5,3} and {5,5,4} were 
each split into two to provide two distinct a_base_elimination lookup tables. 

TableQ]shows the percentage of c base and / coefficients combinations eliminated 
for all c bases and for / coefficients ranging from 2 to 100,000. The elimination per- 
centages could have been calculated based on all possible / coefficients (by ranging the 
coefficient residue from to modulus- 1), but then it would be impossible to exclude 
coefficients that were not z-th power free. The current table is more representative of 
actual searching scenarios, and the ratios between the two calculation methods differ 
by only a few thousandths of a percent. 

B Calculation of / coefficient, c base, and a base Ranges 

This appendix details the logic to calculate the valid range of values for the loops in 
Algorithm Q] in £13.21 Let V = the desired minimum Pegg Value. For exponent sets 
{x, x, z}, the absolute minimum for the three bases is established by £11.6l which states 
that the Pegg Value can be no higher than the smallest original base to the highest 
exponent. For the exponent set {3,4,5}, Lemma |4] necessitates the use of only the 
three coefficient possibilities listed below (and also provides that the resultant equa- 
tion will not be further multiplied by a number to re-associate mm(A, B, C) to a 
base with a higher exponent). For these {3,4,5} scenarios, ill. 61 details how to de- 
termine the absolute minimum for the three bases. We summarize this information 
here: 

if exponent set = {x, x, z}, z > x then 

dminl ^ 1 

else if exponent set = {x, x, z}, z < x then 

bminl V 
Cminl ^ 1 
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else if exponent set = {3,4,5} then 

if coefficient is to be associated with the base to the 3rd power then 

[base to the 3rd power] m j n i <— V/8 

[base to the 4th power] m ;„i <— V/2 

[base to the 5th power] m i„i <— V 
else if coefficient is to be associated with the base to the 4th power then 

[base to the 3rd power] mi „i <— V/4 

[base to the 4th power] mi „i <— V/2 

[base to the 5th power] m ;„i V 
else if coefficient is to be associated with the base to the 5th power then 

[base to the 3rd power] m i„i V/2 

[base to the 4th power] m i„i <— V 

[base to the 5th power] mi „i <— V 
end if 
end if 

The base minimums may be further increased by considering the permutation type and 
combination of other bases. 

B.l Calculating c base Range 

For a given / coefficient, the following logic is used to determine the valid minimum 
and maximum c bases that allow for the resultant equation to fit inside the desired 
range. 

Let ./V = the smallest equation multiplier (based on /) to convert the Original Pegg 
Equation to a resultant equation 

Let Smin = the desired minimum resultant equation size 
Let S m ax = the desired maximum resultant equation size 

B.l.l Permutation Type ax_minus_cz, c base Minimum 

a x - fc z = b y 
fc z +b y = a x 
N(fc z + b y ) = Na x = resultant equation size 

N(fc z mm2 + W) > S mm 

This information does not increase the minimum possible value for c because b can 
always be large enough to make N(fc z nin2 + b y ) > S m i n - The minimum c base is 
Cmini as set initially. 
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B.1.2 Permutation Type ax_minus_cz, c base Maximum 



a x - fc z = b y 
fc z + b y = a x 
N(fc z + b y ) = Na x = resultant equation size 

^(/ c mail + Kninl) — ^max 

f C maxl + bminl — Smax/N 

f c max — Smax/N — 6^jj nl 
c max ^ {Smax/N — & m j n i)/ / 



i/(S max /N-b y mml )/f 



B.1.3 Permutation Type cz_minus_ax, c base Minimum 

For a given /, c must be large enough so that the resulting equation will be > S„ 

fc z -a x = b v 

fc z = a x + b y 
Nfc z = N(a x + b v ) = resultant equation size 

N J ' C-minl — Smin 

c min1 ^ S m .i„ /(Nf) 



Cmin2 — 



</S mm /(Nf) 



and c must be large enough so that fc z > the sum of the minimum a x and b v : 

fc z -a x = b y 

fc z =a x + b y 

fr z > n x A-h v 

J — mini ' mini 



minZ 
Cmin3 



— ( a minl + Kninl)/ f 



\J ( a minl + timinl)/ f 



The minimum c base is the maximum of the above two minimums and the c m i n i as 
set initially. 

Cmin — rn.cix(c m 2nl 3 ^mm2 j Cmin3 ) 

B.1.4 Permutation Type cz_minus_ax, c base Maximum 



fc z ~a x = b v 

fc z = a x + b y 
Nfc z = N(a x + b y ) = resultant equation size 

N f c max — Smax 

c z max < S max /(Nf) 

<JS max /{Nf)\ 
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B.2 Calculating a base Range 

For a given / coefficient and c base, the following logic is used to determine the valid 
minimum and maximum a bases that allow for the resultant equation to fit inside the 
desired range. 

Let N = the smallest equation multiplier (based on /) to convert the Original Pegg 
Equation to a resultant equation 

Let S m in = the desired minimum resultant equation size 
Let S max = the desired maximum resultant equation size 

B.2.1 Permutation Type ax_minus_cz, a base Minimum 

For a given / and c, a must be large enough so that the resulting equation will be 
> f ■ 



fc 



a x = fc z + b y 



Na x = N(fc z + b y ) = resultant equation size 



Nn x > <? 

a min2 — Smin/N 
a min2 = \J Smin / N 



and a must be large enough so that a x > the sum of fc z and the minimum b v : 



- fc z = b y 

a x = fc z + b y 

a min3 — f C + ^rninl 



\J fc Z + &Lnl 



The minimum a base is the maximum of the above two minimums and the a m i n \ as 
set initially. 

B.2.2 Permutation Type ax_minus_cz, a base Maximum 



a x - fc z = W 

a x = fc z + W 
Na x = N(fc z + b y ) = resultant equation size 
Nn x < S 

a max — S m ax/N 



Q"max — ySjnax/N 
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B.2.3 Permutation Type cz_minus_ax, a base Minimum 

For a given / and c, a must be large enough so that the resulting equation will be 

> <? ■ • 

fc z - a x = W 
a x + b y = fc z 
N(a x + b y ) = N(fc z ) = resultant equation size 

N(a x + V) > S mm 

This information does not increase the minimum possible value for a because b can 
always be large enough to make N(a x + b y ) > S m i n . With this permutation type 
a x + b v = fc z . Ifx^y, then this information does not increase the minimum possible 
value for a because b can always be large enough to make a x + b y = fc z . In this case: 

drain ^mml 

But when x = y one of a x and b v must be > fc z /2, and the other must be < fc z /2. 
As established in ^3.21 for searching efficiency we label the larger value a x and the 
smaller value b v . When x = y, this increases the minimum a base to: 

<n„2 > f0 Z /2 



drain — ruax(tt m ^l ; Q'min'l) 

B.2.4 Permutation Type cz_minus_ax, a base Maximum 

fc z -a x =b y 

a x = fc z - b y 

n X < f r z _ h y 
raax — J ° w minl 



\J fc z — b v minl 



B.3 Calculating Valid / coefficients 

For the exponent set {3,4,5}, Lemma [4] details the valid coefficients to be checked. 
This section establishes the valid coefficients for exponent sets {x, x, z}. 

Per i jl.61 the Pegg Value of a resultant equation can be no greater than the smallest 
original base to the highest exponent. We can use this information to determine the 
largest original coefficient that could result in an equation having a minimum desired 
Pegg Value. Given exponent set {x, x, z}, 

Let V = the desired minimum Pegg Value 

Let Smax = the desired maximum resultant equation size 

Let R = the resultant coefficient as part of a resultant base to the highest exponent 
Let H = max(x, z) 

(RraaxV^) — ^max 



RraaxV ^ 



Umax ^ \/ &max /V 
Rmax ySraax jV 
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If R > R m ax, it would limit the original base associated with the highest exponent 
to a value smaller than the desired Pegg Value, thereby precluding the resultant equa- 
tion from having a Pegg Value > V. But our search algorithm does not directly concern 
itself with resultant coefficients. It loops through a range of original coefficients. We 
must answer the question "What is the highest original coefficient that would generate 
an R < R max T 

Consider an Original Pegg Equation {x, x, z} with a single coefficient / associated 
with the base to the z-th power. With {x, x, z}, z < x, the two bases with the highest 
exponent will have a resultant coefficient that is based solely on the multiplier N. 
With {x, x, z}, z > x, the sole base with the highest exponent will have a resultant 
coefficient that is based on the product of / and the multiplier N. 
Define the function cvt(cc, z, v p (f)) for x > 3; z > 3; < v p (f) < z 

if z < x, cvt = smallest^ (-Up (/))/ max(a;, z) 

if z > x, cvt= (smallcst 9 (w p (/)) + v p (f))/ max(x, z) 
using the function: 

smallest g (i;p(/)) = smallest q where q = (mod x),q= —v p (f) (mod z) 
For a p | /, cvt(x, z,v p (f)) returns the power of p as represented in the resultant 
coefficient of the base(s) to the highest power. 

Table 7: Power of p as represented in R (the resultant coefficient of the base(s) to 
the highest power). Output of cvt function for all exponent sets {x, x, z} with each 
exponent < 5. A subset of the information in Table [2] 



Original 
Coefficient 
associated 



Exponent 


with base to 


Highest 




Up(coeff) 


set 


this exponent 


exponent 


1 


2 3 4 


{4,4,3} 


3 


4 


2 


1 1 


{5,5,3} 


3 


5 


1 


2 | 


{3,3,4} 


4 


4 


1 


2 3 


{5,5,4} 


4 


5 


3 


2 1 


{3,3,5} 


5 


5 


2 


1 3 2 


{4,4,5 } 


5 


5 


1 


2 3 4 



With the cvt function to show how a prime dividing an original coefficient is rep- 
resented in the resultant coefficient of the base(s) to the highest exponent, we can 
answer the question "What is the highest original coefficient that would generate an 
R < Rmaxl" The highest possible original coefficient must be < R ma x T where T = 
maximum ratio of v p (f) to cvt(a;, z, v p (f)). 

Consider the exponent set {5,5,4}. For p | / with v p (f) = 3, the entire equation 
must be multiplied by p 5 (as shown in Table [TJ. This is represented as part of the 
resultant coefficient of the bases to the highest power as p 5 ^ z and the cvt function 
shows this with the return value of 1 (representing \og p (p 5 ^ 5 ). So with v p (f ) = 3, an 
original coefficient of p 3 generates an R of p. This 3:1 power ratio is the highest ratio 
for the {5,5,4} exponent set. To produce an R < Rmax, the original coefficient of a 
{5,5,4} equation must be < R m ax 3 - 

For {x, x, z} with all exponents < 5, {5,5,4} has a highest ratio of 3, {4,4,3} 
and {3,3,5} each have a highest ratio of 2, and the remaining exponent sets {5,5,3}, 
{3,3,4}, {4,4,5} each have a highest ratio of 1. 

At this point, the list of possible original coefficients is all integers > 2 and < 
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Rmax ■ The strategy is to create a "valid original coefficient" boolean array of these 
integers, initially marking all the entries as valid. Then we will perform several steps, 
each time possibly ruling out candidates. 



Step 1: Per i jl.21 mark as invalid any coefficients that are not z-th power free. 

Step 2: Mark as invalid any coefficients that have a smallest multiplier N that is 
sufficiently large that it precludes an original equation that has both: 

• a base to the highest exponent that is > the desired Pegg Value 

• that converts to a resultant equation < the maximum equation size 

For example, exponent set {3,3,5} requires a coefficient with v p (f)= 2 to have the 
entire equation multiplied by p 3 , but a coefficient with v p (f)= 1, 3, or 4 requires a 
larger multiplier. Consider the search for Pegg Values > 63742 in {3,3,5} equations 
< 2 88 The highest possible resultant coefficient is 3, as 



Rn 



V / 2 88 /63743 



Let T = 2 = the maximum {3,3,5} ratio of v p (f) to cvt(x, z, v p (f)). The maximum 
original coefficient is given by: 



max _ Umax 
2 



Omax — R 

(1 < 3 

^max _ "J 
^max _ J 

Examining each possible coefficient > 2 and < 9, only the coefficients 4 and 9 can 
product a resultant equation < 2 88 . Step 2 marks as invalid the coefficients 2, 3, 5, 6, 
7,8. 



Step 3: for each coefficient still marked as valid, for the given permutation type, 
calculate the minimum and maximum possible c bases that could result in an equation 
with a desired Pegg Value and within the resultant equation range. If a coefficient 
produces a minimum c base > its maximum c base, the coefficient is marked as invalid. 
For example, consider the search for Pegg Values > 63742 in {5,5,3} equations < 2 100 
of permutation type cz_minus_ax. The highest possible resultant coefficient is 16, as 



Rmax — 

Rmax — 
Rmt 



ySmax /V 



^^21^/63743 



16 



The coefficient 16 is marked as invalid in Step 1 because it is not z-th power free. 



For the original coefficient 15, the minimum c base is given by i jB.1.31 



Crnin3 



\J fenl + & L„l)// 



Cmin ' ^minS 51963742 



v/(63743 5 + 63743 5 )/15 
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For the original coefficient 15, the maximum c base is given by i jB.1.41 



<JS max /{Nf)\ 

^/2 100 /((3 5 * 5 5 ) * (3 * 5))J 
48100619 

So coefficient 15 is excluded as it has no valid c bases. All the other coefficients that 
passed the tests in step 1 and step 2 are valid, as each coefficient has a minimum c base 
< its maximum c base. 

After step 3, all the coefficients still marked as valid are to be checked. 



Crnax — 

Crnax — 
Cmax — 
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